Comment on "Temperature dependence of the Casimir force for lossy bulk media' 
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Recently Yampol'skii et al. [Phys. Rev. A 82, 032511 (2010)] advocated that Lifshitz theory is not 
applicable when the characteristic wavelength of the fluctuating electromagnetic field, responsible 
for the thermal correction to the Casimir force, is larger than the size of the metal test bodies. It was 
claimed that this is the case in experiments which exclude Lifshitz theory combined with the Drude 
model. We calculate the wavelengths of the evanescent waves giving the dominant contribution to 
the thermal correction and we find that they are much smaller than the sizes of the test bodies. 
The opposite conclusion obtained by the authors arose from a confusion between propagating and 
evanescent waves. 



PACS numbers: 31.30.jh, ll.10.Wx, 73.61.At 

It is the subject of a considerable body of literature 
that theoretical predictions for the thermal Casimir force 
between lossy metal plates described by the Drude model, 
based on Lifshitz theory, are in disagreement with experi- 
mental data (see, e.g., review [l[). In Ref. 0] an attempt 
is undertaken to explain this contradiction by arguing 
that Lifshitz theory is indeed inapplicable to test bod- 
ies of finite size, such as those used in the experiments. 
According to Ref. 0, the thermal electromagnetic fluc- 
tuations responsible for the predicted large thermal cor- 
rection excluded by several experiments, have a char- 
acteristic wavelength which is larger than the size of the 
test bodies used in the experiments. On this basis a 
conclusion is made that the predicted correction can be 
observed experimentally only for sufficiently large metal 
bodies. Below we show that the wavelengths of the fluc- 
tuations contributing to the large thermal correction en- 
gendered by the Drude model, are in fact much less than 
the sizes of test bodies used in related experiments. Be- 
cause of this, the purported explanation of the contra- 
diction between experiment and theory in Ref. @ is in 
error. We argue that the considerable overestimate made 
in Ref. 0] of the wavelengths of the contributing fluctu- 
ations, was the result of a confusion between travelling 
(propagating) and evanescent waves. 

The frequencies and wave-vectors of the fluctuating 
electromagnetic field giving a major contribution to the 
thermal correction to the Casimir force can be found 
using Lifshitz formula written in terms of real frequen- 
cies. In modern notation, the thermal correction to the 
Casimir force per unit area, between two parallel semis- 
paces at temperature T separated by a gap of width I, 



can be represented in the form 0, [1[ 
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Here, k± = \k±\ is the magnitude of the projection 
of the wave vector onto the boundary planes, u) is the 
wave frequency, ks is the Boltzmann constant, and 
q 2 = q 2 (uj,k±) = k 2 ^ — lo 2 / c 2 . The reflection coefficients 
for two independent polarizations of the electromagnetic 
field (transverse magnetic, s = TM, and transverse elec- 
tric, s — TE) are given by 
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and e[w) is the frequency-dependent dielectric permit- 
tivity of the material of the semispaces. Equation ([1]) 
coincides with Eq. (3) of Ref. 0], after correcting one 
misprint contained there (in the exponent in the Boltz- 
mann factor on the right-hand side of Eq. (3) the factor 
of two should be erased; an analogous misprint should be 
corrected in Eq. (5) of Ref. @). 

In Ref. [f| it was shown that if the material of the 
semispaces is described by the Drude model 
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where u> p is the plasma frequency and v is the relaxation 
parameter, the major contribution to i^ ra( j is given by TE 
evanescent waves. For example, for Au semispaces with 
hoj p = 9eV and v = 5.32 x 10 13 rad/s at a separation 
I = 162 nm, the TE evanescent waves contribute about 
99.7% of the thermal correction. This contribution can 
be denoted as F r ° a v d an TE . In Eq. Q, the quantity i^ a v | n TE is 
obtained by taking the term with s = TE, for frequencies 
oj varying in the interval from to ck±, for which the 
quantity q is real. Even though the concept of evanescent 
waves is never mentioned in Ref. Q, the contribution 
of TE evanescent waves is actually reproduced by the 
quantityin the first pair of square brackets in Eq. (3) 
of Ref. [2| , integrated over imaginary values of p ranging 
from iO to ioo. The same contribution can be physically 
interpreted in terms of interaction of eddy currents @i[8f. 

As it was also shown in Ref. @, at short separations 
between two semispaces described by Eq. Q, the fre- 
quencies uj giving a dominant contribution to F^ 11 ^ 
satisfy the inequality oj < v(<jj c /uj p ) 2 where ui c — c/(2l) is 
the characteristic frequency. This result was qualitatively 
confirmed in Ref. [2| where the frequencies contribut- 
ing to the quantity F^ n TE — i^^W were found 
to satisfy the inequality oj < v (at I 



100 nm, 



We note that the term F r ° a v d an TE 



it holds 
, which 
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is subtracted from F r ( a ^ I i rE in Ref. Q , represents a negli- 
gibly small thermal effect that results once the material 
for the semispaces is described by the plasma model, and 
it does not influence any of the obtained conclusions. 

It is important to realize that the characteristic wave- 
length of the evanescent waves giving the largest contri- 
bution to the thermal correction is determined, however, 
not by the frequency spectrum of F^ n TE , but rather 
by its wave-vector spectrum. In order to determine the 
latter spectrum numerically, we have recast the quantity 
-^rad "te m the following equivalent form in terms of di- 
mensionless variables 
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where 
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Here, the new variables are defined as 
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We have computed the range v\ < v < vi of the vari- 
able v which contributes 90% of -F'^'xe; at the experi- 
mental separation I = 162 nm. For v = 5.32 x 10 13 rad/s, 
we found v\ — 0.26 and i>2 = 3, while for v — 10 10 rad/s 
we obtained v\ = 0.28 and V2 = 3. Thus, independently 
of the values of the relaxation parameter v, the dimen- 
sionless quantity v contributing to ^^"te ' s always of 
order one, and therefore q is always of order l/l. 

The wave- vector of an evanescent wave is given by the 
expression 
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and its wavelength is determined as 
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Keeping in mind the definition of q, we then obtain that 
for the most contributing wave- vectors it holds k\> q 2 ~ 
l/l 2 , in such a way that the corresponding wavelengths 
satisfy the inequality 



A < 2vrZ 
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In the experiments aiming at measuring the Casimir force 
between a sphere and a plate, these wavelengths are al- 
ways much smaller than the characteristic size L of the 
part of the sphere surface 



Lw2V r R2 



(12) 



which contributes to the force. For example, in the ex- 
periment of Ref. Q the sphere radius is R — 150 /*m, 
and the separation distances vary from I = 162 nm to 
I = 750 nm. For such values of R and I, the inequality 
Amax = 27rZ < L, i.e., I < 2R/-K 2 s» 30/Ltm is satisfied 
with large safety margins, for all the separations consid- 
ered. In fact, the relevant contributing wavelengths are 
smaller than the sizes of the bodies, in all other experi- 
ments measuring the Casimir force performed up to date 
as well [l|. 

The opposite conclusion obtained in Ref. [2| is caused 
by a confusion between propagating and evanescent 
waves. Starting from a qualitatively correct inequal- 
ity for the contributing frequencies to < v, the authors of 
Ref. [2| used the following relation between the frequency 
and the period 
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In terms of these variables the TE reflection coefficient is 
given by 
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to obtain the estimate A > 2-ncjv for the wavelengths 
of the fluctuations contributing to ^Xte- Thereafter, 
it was concluded that Lifshitz theory is only applicable 
if the size of test bodies L ^> 2-ncjv, i.e., v » 2nc/L 
(Eq. (9) in Ref. @)- The problem with this argument, 
though, is that Eq. (TT3l is valid only for travelling (propa- 
gating) waves in vacuum. In this case the two definitions 
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of the wavelength A = 2ttc/lu = 2ir/\k\ coincide. Un- 
fortunately, in the case of evanescent waves, which do 
not propagate and are more similar to standing waves, 
Eq. (fT3")) does not hold, and the wavelength has no rela- 
tion to the frequency. If instead of using Eq. (fT"3]) . the 
authors of Ref. had considered the characteristic val- 
ues of their parameter x = 2ipu)l/c (where in our no- 
tation p = —iqc/uj) to determine the most contributing 
wavelengths, our result A < 2irl would have been repro- 
duced. Indeed, as shown in Ref. 0, x = 2 1 q ~ 1 leading 
to q ~ 1/(2 m qualitative agreement with our esti- 
mate pip. Bearing in mind that for evanescent waves 
the frequency is unrelated to the wavelength^the second 
inequality, co < ksT/h, considered in Ref. Q does not 
lead to any constraint on the size of bodies L. For the 
same reason the results of numerical computations pre- 
sented in Figs. 1 and 2 of Ref. [2[ do not contain any 



information concerning the role of finite sizes of the test 
bodies in calculations of the thermal Casimir forces. 

To conclude, the problem of the disagreement between 
the experimental data of several experiments and the the- 
oretical prediction of the thermal effect in the Casimir 
force, obtained by using Lifshitz theory in combination 
with the Drude model, remains unsolved. 
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